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o : 

, Abstract. Elastic scattering probes directly the interaction potential. For 

weakly interacting condensates this potential is given by the condensate density. 
tJ} We investigate how the differential and total cross sections reflect the density. In 

^ ' particular, we have determined which signatures the Thomas Fermi approximation 

^* leaves in contrast to an exact solution for the condensate wave function within 

the Gross-Pitaevskii theory. 

t> : 

t-H ' After the pioneering experimental realization of Bose-Einstein condensation in 

magnetically trapped atomic gases |[ || these condensates are now routinely 
generated in many laboratories around the world. Consequently the interest is 
shifting from the production and structure of the condensates to their dynamics, i.e., 
their interaction with matter and radiation. One fundamental form of interaction is 
scattering. While it is not entirely clear up to now (see, e.g., how to perform 
experimentally a "standard" scattering experiment (that is to generate projectiles 
with the appropriate extremely low absolute kinetic energies) first theoretical studies 
C"| ■ have explored possible effects || . 

Oh, The purpose of our work is to investigate the signatures of the target distribution 

+-> ' (i.e. the spatial condensate density) in the elastic cross section. In particular, 

we are interested in the differences between an exact description of the density 
in the framework of the Gross-Pitaevskii theory and the additional Thomas Fermi 
approximation which allows one to determine the condensate density analytically in 
k/ the limit of a dense and strongly repulsive atomic gas. 

. ,— i \ We will study these effects in the first Born approximation (see also [|| ) a natural 

starting point for collisionally probing a system at high impact energies (in terms 
of the binding energy of the target, i.e. the condensate). We start from the same 
setup for the weakly coupled many-body Bose system as in and decompose the 
Hamiltonian from the beginning into 

H = H T + Hp + U (1) 

where Ht stands for the target part containing the condensate, Hp represents the 
projectile Hamiltonian of the impinging atom and U describes its interaction with the 
target atoms. The target Hamiltonian reads 

Ht= J d 3 r&(r) (~^A r + V ex + V a (r)\ *(r) , (2) 

where fy(r) is the field operator for destruction of bosons at r, Vex = mu> 2 r 2 /2 is the 
harmonic trapping potential. The contact interaction between the atoms of mass m 
is expressed through the s-wave scattering length a s as 

V a (r) = lj d 3 r'¥ (r')V(r, r')*(r') (3) 
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with 

V(r,r') = 47i7i 2 a s /m<5(r™r'). (4) 

For simplicity we consider projectile atoms of the same sort as the condensate atoms. 
Hence, the interaction potential 

U(r p ) = 2V a (r p ) (5) 

of the condensed atoms with the projectile atom is of the same form as the potential 
V a (r) within the condensate, given in (^), where r p is the position of the projectile 
atom. The factor two has its origin in the single particle nature of the interaction of 
the scattering atom with the condensed atoms. 

In the weak coupling limit we may write the field operator as [0 

*(r) = iV 1/2 * (r) + 5* (r) (6) 

where "Jo is a c-number whose squared modulus represents the normalized density 
of the condensate at position r. Minimizing the grand canonical Hamiltonian 
Kt = Ht — [iNt, where Nt is the number operator of the trapped atoms, with 
respect to the functional [^o] to zeroth order in 6^(r) gives the time-independent 
Gross-Pitaevskii equation for the condensate wave function "to 

-^A r + Vex + N ^h 2 a s /m\m Q \ 2 - fj) * = 0. (7) 

The eigenvalue /j, is the chemical potential. With (Q) and 'J'o from (Q) one obtains the 
elementary quasiparticle excitation modes and energies by diagonalizing Kj> to second 
order in S^f. This is the "standard" Bogolubov approximation (?], 

However, in the present context we are only interested in elastic collisions. In 
leading order of the number of condensed atoms No the excitations do not contribute 
to these processes as we will see shortly. 

The asymptotic scattering states are eigenstates of Ht and Hp respectively, i.e. 
condensate eigenstates and |<I>/} as described above, and plane waves for the 
projectile atom with asymptotic initial and final momenta and k/, respectively. 
With equation (Q) the scattering amplitude in first Born approximation reads then 

F(M) = / <P™~ ik, '(*lMr)\*i)e <k " (8) 

with the cross section defined as 

Inserting (0) into (@) and keeping only leading terms in iVo one obtains the elastic 
scattering amplitude 



F e\( 



) = -2a s J d 3 re- 4qr 7V |* (r)| 2 (10) 

where q = ky — k^ is the vector of momentum transfer and no energy has been 
exchanged between initial and final state, i.e. <&i = $/ and |kj| = |k/| = k. 

From (jl^) one can see that the elastic cross section probes directly the condensate 
density l^ol 2 - For structural properties of repulsive atomic gases with large particle 
number iVo it is often sufficient to determine m the Thomas-Fermi approximation 
[§ as applied, e.g., in @. 
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Its justification is most easily seen if one scales the unit length and energy 
in (^) such that the only remaining parameter appears as an effective mass T. 
Before doing so, we separate the relevant radial s-wave part by the substitution 
^o( r ) = M(r)/r(47r) -1 / 2 and obtain in energy and length units of the trapping 
potential, namely fkv and a u — [fi/lmu)] 1 ' 2 , 

where 

T = N a s /a u . (12) 

Scaling the length in ( |li"| ) according to r = r 1 / 4 ^, dividing ([ll]) by T 1 / 2 , and defining 
u(x) — u(r) leads to 

1 d 2 x 2 u{xf 



2Tdx 2 ' 2 ' x 2 ^) = > ( 13 ) 

where ft = /i/huT 1 / 2 . For large T the kinetic energy becomes small and is neglected 
in the Thomas Fermi approximation (TFA). This leads to the solution 

u(r)/r = ( M /r(l - {r/Rf)f 2 , (14) 

where n(T) = fia;/2(15r) 2 / 5 is obtained from the normalization condition J u 2 (r)dr = 
1. It defines simultaneously the cutoff R(T) — a w (15r) 1 / 5 which obeys the relation 

muj 2 „n 

fl =—R 2 . (15) 



The solution of (|llj) must be obtained numerically, e.g., by imaginary time propagation 
with the Split Operator Method 0. 

From the construction it is clear that the TFA is only applicable for large and 



positive r which implies according to (12) a positive scattering length a s > and 



under realistic experimental conditions a large particle number No 3> 1. This is 
illustrated in figure [l] where the (normalized) solution ^o(r) of (Q) is compared to 
the TF-approximate solution ([l4|). Significant differences appear for T < 10, for the 
order parameter (part a) as well as for the chemical potential (part b). However, 
from problems involving a simple atom as a target, one knows that elastic scattering 
at higher impact energies probes the target distribution to some detail. Hence, it 
is a priori not clear to what extent the TFA describes cross sections well, even in 
a parameter regime for T 3> 1 where it works well for structural properties. For 
total elastic cross sections the best agreement, as expected, is obtained at low impact 
momentum k (figure |J). The functional dependence on k is the same for the exact 
solution and the TFA approximation for large k and can be derived from the TFA 
approximation to be a(kR 3> 1) oc (Rk)~ 2 , where R = a w (15r) 1 / 5 is the condensate 
radius in the TFA approximation, see (jl4|). However, the exact solution and the TFA 
differ by a factor which grows with decreasing T (figure ^a) . 

On the other hand, the total cross section for the TFA in the Born approximation 
can be mapped to a single universal curve if one scales a = cr/T 2 and expresses the 
momentum in the dimcnsionlcss varibale k — kR where R is the cutoff parameter of 
(fLq). The numerical solution does not have a well defined cutoff parameter (see figure 



hja) . However, we may use (15) to define the cutoff R^ from the numerically obtained 
chemical potential [i. With this scaling the cross sections from the numerical order 
parameter fall onto the same universal curve as the TF cross sections (||d) . 
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Considerably more disagreement between the TFA and the exact solution is seen 
in the differential cross section shown as a function of momentum transfer in figure [|. 
One clearly recognizes the exponential versus the algebraic decrease in the absolute 
magnitude and the different oscillatory behavior. Both effects can be understood 
analytically. The "soft edge" of the numerically obtained condensate distribution 
leads to an exponential decrease of the cross section with increasing q. A similar 
effect occurs in the photoionization and the Born cross section for scattering of metal 
(jellium) clusters due to a soft edge in the electron distribution [jlT], |l2). On the other 
hand, the "hard edge" of the analytical TFA distribution (|lj) produces via Fourier 
transform the (unphysical) algebraic decrease of the differential cross section. (In the 
photoionization of clusters this would correspond to a box-like electron distribution 
with a sharp cut-off, Jll[). For the oscillatory behavior we note that the exact cross 
section oscillates asymptotically for large q and T with the same frequency as the TF 
approximated cross section given by the condensate radius R. 

Overall we conclude that the TFA offers also for scattering problems a reasonably 
good approximation as long as T is not too small. Significant deviations from the exact 
solution occur only for large momentum transfer q for which the cross sections are so 
small that they will not be experimentally accessible. 
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Figure 1. Properties of the numerical {solid, ([ll])) and Thomas Fermi {dashed, 
(114)) solution of the Gross-Pitaevskii equation. Part (a) presents the order 
parameter ipo{r) normalized to unity at r = for Y = 1/10, 1, 10, 100, 1000 (from 
left to right). Part (b) shows the chemical potential as a function of T. 
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Figure 2. Total elastic cross sections as a function of V for k = 0.2, 1.2, 5 from 
top to botton. Coding of the lines as in figure hi 
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Figure 3. (a) Total elastic cross sections for the same V values as in figure |lja 
from top to bottom. Part (b) shows the cross sections of (a) in scaled coordinates, 
see text. The symbols (+,X,o,D,0) correspond to Y = 1/10, 1, 10, 100, 1000 and 
the dashed line is the universal TF curve. 
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Figure 4. Differential elastic cross sections as a function of momentum transfer 
q for r = 1/10, 10, 1000 from bottom to top in both the group of numerical cross 
sections (solid) and analytical cross sections (dashed). 



